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Abstract 



^ c"| In this paper, we first introduce multifrational Riemann-Liouville Brownian sheets. 



Then, we show a result of approximation in law of the multifractional Riemann- 
Liouville Brownian sheet. The construction of these approximations is based on a 
sequence of I.I.D random variables. 
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CN . 1. Introduction 

CN ' A one-dimensional fractional Brownian motion (FBM) = t € K+} with Hurst 

index H G (0, 1) is a real-valued, centered Gaussian process with covariance function given 
by 



E 



t H {t)t H {s) = \\\s\ 2H + \t\ 2H - \t-s\ 2H ], s, t £R+. (1.1) 



2 

It was introduced, as a moving-average Gaussian process, by Mandelbrot and Van Ness 



Fractional Brownian motion possesses interesting properties such as self-similarity of 
order H € (0, 1), stationary increments and long-range dependence (when H > i), which 
make it a good candidate for modeling different phenomena in, for example, finance and 
telecommunication. However, it does not represent a casual time-invariant system as 
there is no well-defined impulse response function. Barnes and Allan [2] introduced the 
fractional Riemann-Liouville (RL) Brownian motion (RL-FBM) Vh = {Vff(i); t G K + } 
based on the following Riemann-Liouville fractional integral, 

Vh ^ = WW~U I (* - u)+~ h dB(u), (1.2) 

1 [H + o j 7R+ 
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where B(u) is a standard Brownian motion, x+ = max{x, 0}, and V is the gamma function. 
Vh represents a linear system driven by white noise with the impulse response function 



* \, . Vh shares with £h many properties which include self-similarity, regularity of 

sample paths, etc.- with one notable exception that its increment process is non stationary. 
But it was shown by Lim [14] that the increment process satisfies some weaker forms of 
stationary. This provides some flexibility necessary for practical applications (see Lim 
|14j for more information). However, this model may be restrictive, due to the fact that 
all of its regularity and fractal properties are governed by the single Hurst parameter H. 
In order to overcome this limitation, Lim and Muniandy [16] extended RL-FBM to RL- 
multifractional Brownian motion (RL-MBM). Roughly speaking, RL-MBM is governed by 
a Hurst function H (i) with certain regularity in place of the constant Hurst parameter 
H in RL-FBM. We refer to Lim p3] and Lim and Muniandy \15\ [T6] for properties of 
RL-MBMs. 

There are two typical multiparameter extensions of FBMs. One is the Levy's fractional 
Brownian random field with parameter H £ (0, 1) (see Ciesielski and Kamont [6]) : a 
centered Gaussian process Y with covariance function given by 



i 



E[Y(t)Y(s)] = \ 



\s\\ 2H + \\t\\ 2H - \\t - s\\ 2H 



, s,te m+, 



where || • || denotes the Euclidean norm. 

The other multiparameter extension of the fractional Brownian motion is the anisotropic 
fractional Brownian sheet introduced by Kamont [11]. Since then, fractional Brownian 
sheets have been studied extensively as a representative of anisotropic Gaussian random 
fields. See, for example, Wu and Xiao [18] and the references therein for further infor- 
mation. Still, the regularity of fractional Brownian sheets does not evolve in the d— 
dimensional time parameter t G Ml. 

To model the anisotropic Gaussian random fields whose regularity evolves in time, such 
as images, Ayache and Leger [1] and Herbin [9] introduced so-called multifractional Brow- 
nian sheets (MFBS) in terms of their moving average representations and harmonizable 
representations, where the constant Hurst vector of fractional Brownian sheets is substi- 
tuted by Hurst functionals. We refer to Ayahce and Leger p], Herbin [9], Meerschaert, 
Wu and Xiao [13] and the references therein for more information. 

Now, we are ready to introduce the multifractional Riemann-Liouville Brownian sheet 
(RL-MBS). 

Definition 1.1 Let H{t) = (Hi(t), ■ ■ ■ ,Hd(t)) be a measurable function in t € with 
values in (0, l) d . A real-valued RL-multifractional Brownian sheet X = {X(t); t € 
M^} with functional Hurst index H(t) is defined as the following Riemann-Liouville (RL) 
fractional integral 

X(t)= [ Wti-u^f^Widu), (1.3) 

M t=i 

where W{u) = {W(u); u £ R+} is a (i-parameter Wiener process. 

Some authors have studied weak convergence to multifractional Brownian motions. Dai 
and Li [7] presented a weak limit theorem for the multifractioan Brownian motion based 
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on a Poisson process. By using Donsker's theorem, Dai [8] showed an approximation of 
the RL-multifractional Brownian motion in Besov spaces. On the other hand, Weak limit 
theorems for fractional Brownian sheets have attracted significant interest in recent years. 
For example, Bardina et al. [5] gave a result of approximation in law, in the space of the 
continuous functions on [0, l] 2 , of two-parameter fractional Brownian sheet. Bardina and 
Florit [3] showed a result of approximation in law of the d-parameter fractional Brownian 
sheet in the space of continuous functions on [0, l] d . 

Let us recall some known facts. Let {Z^; k = (fei, &2j ■ ■ ■ , fe<z) €= be an independent 
family of centered identically distributed random variables with variance 1. Wichura [17] 
generalized Donsker's theorem and proved that the process 



where I[ k -i,k)(u • n) = Ifa-i, fci)x-x[fc d -i, k d )( u i n , • • • , u dn), converges in law to a d- 
parameter Wiener process. 

Inspired by Wichura [T7] and Bardina and Florit [3], it is natural to try to approximate 
X in law by 



• For every % G {1, • • ■ , d}, < on < Hi(t) < ft < 1 for all t G M^_. 

• H(t) is 7-Holder continuous, i.e., 

\\H(t) -H(s) || <K\\t-s\\^, Me«t 
where K > is a constant. 

Furthermore, we assume that E|Z^| m < oo for any k € N d and m € N. 

In this paper, we consider [0, l] d C Rl. with the usual partial order. We will prove that 
the sequence of laws in C([0, l] d ) of the processes {X n (t);t G [0, l] d } defined by (|1.5p 
converges weakly to the law of {X(t); t G [0, l] d } given by f)1.3|) . 

Most of the estimates of this paper contain unspecified constants. An unspecified posi- 
tive and finite constant will be denoted by K, which may not be the same in each occur- 
rence. Sometimes we shall emphasize the dependence of these constants upon parameters. 

We end this section with some properties of RL-multifractional Brownian sheets. 

Lemma 1.1 For any s < t G [0, l] d , there exits a constant K > such that 




(1.4) 




(1.5) 



where n (u) = nz Y,km d z k I [k~i, k)( u ■ n). 

In the rest of this paper, we assume that H(t) = (H±(t), • • • , Hd(t)) satisfies the following 
conditions: 



E X(t) - X(s) < K\\t - s 



2H 



where H = min i6 { lj ... )d }{at,7}. 
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Proof: First, we have 



E 



X(t)X(s) 



)} = nil [ / (u- Ul ) H / t] 5( Si - u,)+ (s) 1 

J L ^[o, 1] 



dui 



(1.6) 



By (jl.6p . we can consider the process X = {X(t)} defined from d independent RL- 
multifractional Brownian motions Y^> with parameter Hi(ti) = Hi(t) by 



Y(t)=Y^(t 1 )---Y^(t d ). 



[1.7) 



It is easy to see that it has the same covariance function as the RL-mutifractional 
Brownian sheet X. 
From above arguments, we get 



E 



X(t)-X(s) 



(uf =1 Y^(U) - YW( Si )IlUY®(tij) +■■■ 
+ (n£iY»( ai )yM(t d ) _ uf =1 Y^(s t )] 

Since YW, i = 1, • • ■ ,d are mutually independent and E[FW] 2 < K, we get that 
E[x(t) - X(s)] 2 < Kj2^[Y {i) (U) - Y^(si)} 2 . 
It follows from Dai [5] that 



(1.8) 



i=l 



(1.9) 



E 



yW(ti) - Y^( Si )\ < K\U - Si \ 2H ^ + K\Hi(t) - H^ 2 . 



Since t,s€ [0, l] c 



On the other hand, we have 



\U - sA < 1. 



\Hi(t) - Hi(s)\ < 1. 



By (??) to (|1.12|) . one can easily get that the lemma holds. 



(1.10) 



(1.11) 



(1.12) 



□ 



It follows from the Kolmogorov continuity theorem (see Theorem 2.3.1 in Khoshnevisan 
TUP and Lemma 1 1.1 1 that 



Corollary 1.1 {X(t)} has a continuous modification. 
The proof of Corollary 11.11 is classical, so we omit the proof. 



4 



Approximation to RL-MBS 



5 



2. Main result 

The main result of this paper is the following theorem. 

Theorem 2.1 The sequence of laws in C([0, l] d ) of the processes {X n (t);t G [0, l] d } 
defined by \1.5\) converges weakly to the law of {X(t); t G [0, l] d } given by tl.3\) . 

In order to prove Theorem [27T1 we have to check that the family of laws of the processes 
{X n } is tight. First, we need a technical lemma. 

Lemma 2.1 For any even m G N, there exists a constant K > such that for any 
nonnegative functions fa G L 2 ([0, 1]), i = 1, ■ ■ ■ ,d, 



E 



\fi{ui)\e n {u)du <k] \ fi( Ui )d 

[0, i] d L ~I J L -'" 



i=l 



11 ; 



(2.1) 



Proof: We have 
d 

[Y[fi(ui)]9 n (u)di 



E 



i=l 



■ m 




= E 


[/ 




. J[0, l] d ™ 



[0, l] dm 



m d m 

[ n n /<w)] [ it ^(^a* 1 

3=1 *=1 3=1 
m d m 

[ n n e [ n w)] dui --- ^ m ( 2 - 2 ) 



j=i i=i 



3=1 



On the other hand, we have 



m 
3=1 



dm 

n 2 E 



dm 

n 2 E 



[II ( Yl Z k l \k-\,k)W ■ n) 
3=1 fceN d 

^ Z fc i • • • ^fcmJ[ fe i_ 1) k i)(u l ■ n) ■ ■ ■ /[fem-i, k m )(u m ■ n) 



dm r — / 1 \ / 77 j 

< J^n 2 ^ ^jfci-i,*; 1 )^ •«)■■■ -'[few-l, fc™)^ •«), 



(2.3) 



where 



,4 m = | (A; 1 , • • • ,k m ) G N dm ; for all i G {1, • • • ,m}, 

there exists j G {1, • • • , m} \ {ijsuch that k % = 

since Zj., k G N d are independent and E|Zfc| m < oo for all m G N and k G N d . 
We have 



h^-i, fc 1 )^ 1 ' n ) ' ' ' J [fc m -i, k^)(u m ■ n) < I Dm {u l , 

(fc 1 ,- ,fc m )GA m 



,1i 



(2.4) 



where 



L> m = {(u 1 ,--- ,u m ) G [0, l) dm ; for all i G {!,-•• ,m}, 



there exists j G {1, • • • , m} \ {i} such that for all k = 1, • • • , d, \u\ — uA < 



n ) 
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In fact, we can bound Ijj m by a finite sum of products of indictators. Moreover all the dm 
variables, u\, ■ ■ ■ , it™, appear in each product of indicators, but each variable appears only 
in each one of indicators of each product. Then we can bound (12. 2D by a sum of products 
of the following two kinds of terms: 

• For any i ^ j, i,j £ {1, • • • ,m}, 

„ d d 

h = Kn d / ' [ [f k (ul)f k (u{)} J] J [0 1)04 - uHWduK (2.5) 
J\o,i] 2d , , n 



k=l 



k=l 



For any i=£ j ^ r, i, j,r 6 {1, • • • ,m}, 

d 



Io = Kn 2 



[ [/fcK)/fcK)/K)] 



fc=i 



[] [j [0) i)(K - 4l) J [o, 1)(I4 - «fcl)V ±)(K " 41)] d«WW. (2.6) 

k=l 



In order to prove the lemma, we only need to show 

d 



h<K [ fl[fi(u k )]du 



(2.7) 



and 



h < K 



\[fl{u k )}du 



(2.8) 



First, we show that (I2.7P holds. I\ can be bounded by 



Kn° 



[[/ fc K)] 2 n 7 [o^)(K-4l)d«w 



fc=i 



fc=i 



r. d d 

+Kn d / n [ak)] 2 n / [0 A) (K - 4Dd«w. 

JfO, ll M 1 ™ ; 



fc=l 



fe=l 



(2.9) 



where we have used the elementary inequality 2ab < a 2 + b 2 . 

The two terms in (I2.9D can be done with the same method. So we only deal with the 
first one. In fact, 



Kn d [ n^K)] 2 llVi' 

J \. U ' l i k = l k=l 



(K-<|)d«w 



< k 



[0- fc=l 



II [AK)] 2 * J . 



(2.10) 



From (123]) and ([2~TU1) . we can get that ([221) holds - 
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Now, we prove f|2.8jl . I2 equals 

d 

'[0, l]« L fc=i 



Knf [ [n[/ fe (4)/ [0)i) (i4-4i)]' 

■no. i] 2d L fc =1 n 

/ II [AK)AK)^ [0 , i)(K - i)(K - «JI)]** r 



du l du 3 . 



By the Holder inequality, 



1 1 



h < Kn^Gf x G 2 , 



(2.11) 



where 



V i] 2d ^ L 1 " J 



Ufc-Ufcl 



du l du J , 



and 
G 2 



[0, i] 2 



[0. 



d 

II [/feK)AK) J [o, i)(K - i)(K - 41) 

1 l d fc=l 



dvU'd'u? . 



Using the same method as the proof of (j2.7|) . we can get 



Gt < Kn- d / T /|(u fc )di 
Now, we deal with G2. G2 equals 



(2.12) 



/ n [/fc(<)/fc(«fc)/fc( u fc) j [o, - ki^o, i)(k 2 

•/[0, i] 4d r^; L ™ n 

^[o,i)(l«fc -4lK [0 ,i)(K 2 -<|)]^W 2 duW. (2.13) 

So, G2 can be bounded by 

Kn~ d [ J] [/fcKJAK 1 )/*^*)^, i)(l u fc " 4lK[o, 
-'[0 i] 3d - - 



< 2 



fc=i 



<i^n- d (G 3 + G 4 ), 



du ri du r2 du j 
(2.14) 



where 



7rn 113d L 71 



du ri du r2 du j , 



k=l 
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and 



Ga 



[0, 1]3 



n[^(«j)/*(«3?Kip 1 i)(i«* 



41) 



du ri du r2 du j . 



k=l 



G3 can be bounded by 



k=l 



■Ao, 11 M L ,_, 

d 



du ri du J 



< Kn 



Kn 



[0, 1] 5 



fc=i 



/c=i 



(2.15) 



Using the same method as above, we can get that G4 can be bounded by 



Kn 



n [ft 



(2.16) 



It follows from (12. lip to (|2.16j) that the inequality (12. 8j) holds. This completes the proof. 

□ 

Lemma 2.2 For any even m £ N, there exists a constant K > such that for any 
s,t € [0, l] d with s <t, 



E 



X n (t)-X n (s) 



< K\\t - s 



where H = min^^... ^{a;, 7}. 
Proof: It follows from Lemma |2. II that 



E 



X n (t)-X n (s) 



< K 

< K\\t - s 



[0, i] d 



l((t l -u,)f^-(s t -u 1 )f^Yd l 



i=l 
mH 



□ 



On the other hand, one can easily prove that 
Corollary 2.1 The process X n is continuous for every t 6 [0, l] d a.s. 
Next, we will prove the main result of this paper. 

Proof of Theorem \2.t\ By Bickel and Wichura [5] , in order to prove Theorem 12. 1\ we 
only need to prove that the family of laws of processes {X n (t)} is tight, and identify the 
limit law of any convergent subsequence as the law of the d-parameter RL-multifractional 
Brownian sheet X. The tightness comes from Lemma 12.21 and Bickel and Wichura [5]. 
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Now, we proceed to identify the limit law. It is sufficient to show that for any q G 
N,ai,-- - ,a q GMandt 1 ,--- , t q G [0, l] d , 



E 



i r 

exp{^^ ai X n (t j )} — > E expktJ2 a i X ( tj )} 



j=i 



(2.17) 



as n — > oo. 

For any j G {1, • • ■ , g} and i E {1, • • • , we have 



and 



(^-«i)r ' '€^([o, id, 



(2.18) 



(2.19) 



Therefore, for any j G {1, • • • , g} and i G {1, • • • , d}, we can find a sequence {p? k (ui)}keN 
of simple functions such that 



[o, i] 



and 



.I\2 

' d«j — >• 0, as /c — > oo, 



$ - Ui) 1 ^"' 1 5 - ^(u*) > 0, for all k G N. 
In order to simplify the notation, we define 



(2.20) 



(2.21) 



x 3 - 



[Hp> ik ( Ui )]0 n (u)du, 

[o,il d i=l 



and 



Note that 



E 



' i r 

expjze^a^n^)} -E exp a, 



where 



h 



3=1 



E 



< /x + 7 2 + ig, (2.22) 



E 



exp {i£ " exp ji£ a i X1 n k ] 

3=1 3=1 

q q 
exp [it £ djX(f)) - exp [it £ a 3 X 3 > k } 



3=1 



3=1 
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and 



E 



exp {i£ a 3 X n k } ~ ex P Yl a i Xj,k } 



3=1 

We first deal with I\. By the mean value theorem, 



h < K max <^ E 

i<j<q 



On the other hand, for j € {1, 2, • • • , q} 
E 



x n (ti)-xk k \) 



(2.23) 



E 

'[o, i]° 



(tj-ui) 



j=2 



+ 



■ / [o.i] d L f=2 r=2 J 



< E 



+E 



(*i-«i)? (t ' )-l -pi !fe K)][n^ 

J i=2 
d 



,Hi(t j )-h 



9 n (u)du 



'[o, i] 

On the other hand 



i=2 



i=2 



(2.24) 



E 



[o, i] d 

< # 

< A' 



,(«i)][n(* 



i=2 



(*i-«i)? ( * H -ri lfc 

i — ^ H -^ fc («i))n(*? 



[0, 1]" 



[0, 1] 



i=2 



da 



(^-m)+ l(t) 2 - pj fc («i)) dux 



(2.25) 



where we have used the Holder inequality and Lemma l2.1i 
It follows from flg^ZOJ) and (|2T2l)j) that when fe ->■ oo, 



E 



[o, i] d 



i=2 



0, (2.26) 



where — > denotes uniform convergence with respect to n. 
Repeating the above procedure, we can get 



E 



x n (v)-xt 



— > 0, as k — > oo. 
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Therefore, 



h 0, as k -4 oo. (2.27) 



Similarly, by the mean value theorem, we can get 

h < K max {e|X(^') - X j ' k \\ (2.28) 
Applying the Holder inequality and the properties of the stochastic integral, we can get 



E\X(t J ) -Xi' k \ < K 



J\0, l] d V M M y 



Using the same method as the proof of (|2.27p . we can get that 



(2.29) 



h 0, as k -4 oo. (2.30) 
Due to the result established by Wichura |17j . the law of the process 

r-tl ft d 



/ • • / <?n< 

io io 



u)du 



converges weakly to the law of a d-parameter Wiener process. Then, the linear combina- 
tions of increments of 

rh rtd 

6 n (u)du 

will converge in law to the same combinations of increments of the (i-parameter Wiener 
process. So 

I 3 -> 0, (2.31) 

asn-> oo. 

Combining (f2"T7D . (l2T22l . (12371) . (p3D1) and we can get that Theorem O holds. 

□ 
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